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Abstract 

The perturbative framework is developed for the calculation of the 7r + 7r~ atom 
characteristics (energy level shift and lifetime) on the basis of the field-theoretical 
Bethe-Salpeter approach. A closed expression for the first-order correction to the 
ir + ir~ atom lifetime has been obtained. 
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Recently a number of experiments on the study of hadronic atoms have been 
carried out. Namely, the first estimate of the 7r + 7r - atom lifetime was given in 
Ref. [1]. The measurement of the characteristics of the pionic hydrogen [2] and 
pionic deuterium [3] have been performed. At present, the DIRAC collaborati¬ 
on is preparing the experiment at CERN on the high precision measurement of 
the lifetime of 7r + 7r“ atoms. This experiment might provide a decisive improve¬ 
ment in the direct determination of the difference of the S'-wave tttt scattering 
lengths a[j — ag and thus serve as a valuable test for the predictions of Chiral 
Perturbation Theory [4]. Note that the analogous experiments on the observa¬ 
tion of the ttK, pK atoms are also planned. In the view of these experiments 
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there arises a need in the theoretical framework which would enable one to 
calculate the characteristics of such atoms with a high precision based on the 
ideas of standard model. 

The history of the theoretical study of hadronic atoms has started from 
Ref. [5]. In this paper the relations of the energy level displacement of the 
pn atom due to strong interactions and its lifetime with the strong nN scat¬ 
tering lengths have been established in the framework of the nonrclativistic 
scattering theory. In the following, this approach has been generalized to the 
case of the atom [6,7]. In particular the expression for the lifetime T\ of 

the 7r + 7r~ atom in the ground state was found to be (see also Ref. [5]) 
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where the isotopic invariance of strong interactions was assumed. Here we 
define Am, is the m nr ± — m n o mass difference, and i (0) = 0o = (m^.a 3 /87r) 1 / 2 
is the Coulombic wave function (w.f.) of the pionium at the origin. 


The approach to the study of the problem of hadronic atoms, developed in 
Ref. [5], makes use of the general characteristic feature of the hadronic atoms - 
the factorization of strong and electromagnetic interactions. Namely, since the 
Bohr radii of the atoms composed from the light mesons and nucleons is of 
order of a few hundreds fm, their energy spectrum is almost completely deter¬ 
mined by the static Coulombic potential acting between the constituents. On 
the other hand, the decays of hadronic atoms are governed by the strong inter¬ 
actions which, e.g. for the case of the pionium are responsible for the transition 
of the 7T + 7r _ -pair into 7r°7r°. In the following we shall consider this particular 
case. The formula (1) demonstrates this factorization property explicitly, ex¬ 
pressing the atom lifetime as a product of two factors - the Coulombic w.f. at 
the origin and the strong interaction part, completely concentrated in the 7T7 t 
strong scattering lengths. 


The problem of evaluation of the electromagnetic and strong corrections to the 
basic formula (1) within different approaches has been addressed in Refs. [8]- 
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[14]. For a brief review see Ref. [13]. In this paper within the Bethe-Salpeter 
(BS) approach we have derived the relativistic analogue of the formula (1) 
taking into account the strong interaction corrections in the first order. These 
corrections were found to be of the relative order 10 -3 . It should be stressed 
that the field-theoretical approaches [11,13,14] to the problem, unlike the po¬ 
tential treatment [8,12], do not refer to the concept of the phenomenological 
strong interaction mr potential, which is a source of an additional ambiguity 
in the calculations of hadronic atom characteristics. In the former approaches 
these characteristics are expressed directly in terms of the underlying strong 
interaction (chiral) Lagrangian, and the results can be compared to the expe¬ 
riment, providing the consistent test of the predictions of chiral theory. 


In the present work we suggest a relativistic perturbative framework for the 
calculation of the bound-state characteristics of hadronic atoms (energy levels 
and lifetime). Our framework, based on the BS approach to the bound-state 
problem, is quite similar to the one used in the treatment of the positron- 
ium problem (see, e.g [15]), and most of the methods which are used in the 
latter case can be applied to the hadronic atoms as well. The main purpose 
of this work is to demonstrate the possibility (not only in the potential scat¬ 
tering theory, but in the BS treatment as well) of the clear-cut factorization 
of strong and electromagnetic interactions in the observable characteristics of 
hadronic atoms, avoiding the double-counting problem in the calculation of 
these quantities. One should note that the suggested approach allows to cal¬ 
culate strong and electromagnetic corrections in all orders of the perturbation 
theory. At the present stage we apply the general formalism to the calculation 
of the first-order strong and electromagnetic corrections to the pionium lifeti¬ 
me. The results for strong corrections obtained in Ref. [13] are reproduced in 
these calculations. 


Now we pass to the description of the perturbative framework. The perturba¬ 
tion expansion is performed around the solution of the BS equation with the 
purely Coulombic kernel similar to that introduced in Ref. [16] 
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( 2 ) 


Vb(p,q) = \^(^)iV / ^) 

Here m, denotes the mass of the charged 7r-meson, and w(p) = \jm 2 + p 2 . 
The factor introduced in the definition of the instantaneous Cou- 

lombic kernel (2) enables one to reduce the BS equation with such kernel to 
the exactly solvable Schrodinger equation with the Coulombic potential. Then, 
the exact solution of the BS equation with this kernel is written in the form 

^c(p) = iG 0 (M*-,p) 4y/in(p) 4?r ^^ 0 , ^ c ( p ) = -0 c ( p ) (3) 

v p- + p~ 

where p = m^oi /2 and M* 2 = m 2 (4 — a 2 ) is the eigenvalue corresponding to 
the unperturbed ground-state solution. Go denotes the free Green’s function 
of the 7T + 7r _ -pair. The BS w.f. (3) is normalized in the usual way [17] 

< ipc\N(M*)\ip c >— 1, JV(M*) = (4) 

The exact Green’s function corresponding to the Coulombic kernel (2) can be 
constructed according to Ref. [18] 

Gc(P";p, 9 ) = (2jt) 4 ,5 <4) (p - «)G„(P*; p) + G„(P*; p)T c (£*; p. q)G„(P*; q) ( 5 ) 
where Tq is given by 


T C {E *; p, q) = 1 y/w(p)w(q) 


+ 


vdpp 


p-q ) 2 J D(p; p,q)J 


( 6 ) 


D(p\p,q) = {p-q) 2 p 


4 E*\ m n J V m' y 


where v = a^Jm n / (—4 E*) and E* = ( P* 2 — 4m 2 )/(4m 7r ). 

The full BS equation for the 7r + 7r~ atom w.f. x(p) is written as 


r d^k 

Go\P;p)x(p) = J j^V(P-,p,q) X (q) 


(7) 
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where V (P; p, q) denotes the full BS kernel which is constructed from the un¬ 
derlying (effective) Lagrangian according to the general rules and includes all 
strong and electromagnetic two-charged-pion irreducible diagrams. In parti¬ 
cular, it contains the diagrams with two neutral pions in the intermediate 
state which determine the decay the 7r + 7r - atom into 7r°7r°. Note that in addi¬ 
tion V (P; p, q) contains the charged pion self-energy diagrams attached to the 
outgoing pionic legs (with the relative momentum q), which are two-charged- 
pion reducible. These diagrams arise in the definition of the kernel V (P; p, q) 
because the free two-particle Green’s function is used in the l.h.s. of Eq.(7) 
instead of the dressed one. The c.m. momentum squared P 2 of the atom 
has the complex value, corresponding to the fact that the atom is an un¬ 
stable system. According to the conventional parametrization, we can write 
P 2 = M 2 = M 2 — iMT where M denotes the ’’mass” of the atom, and T is 
the atom decay width. 

Further, we introduce the four-point Green’s function for the 7r + 7r“ —> 7r + 7r~ 
transition which, by definition, obeys the inhomogenuos BS equation 

G(P) = G 0 (P) + G 0 (P)V(P)G(P) (8) 

This function has a pole in the complex P 2 plane at the bound-state energy. 
The relation between the exact solution x(p) and the Coulombic w.f. is 
given by [13] 

<X \=C <^c\ Gc\P*)G(P), P* 2 M* 2 , P 2 M 2 (9) 

where C is the normalization constant. In what follows we assume that this 
limiting procedure is performed with the use of the following prescription [13] 
P* 2 = M* 2 + A, P 2 = M 2 + A, A —■» 0. The validity of Eq.(9) can be trivially 
checked, extracting the bound-state pole in G(P) and using the BS equation 
for i[}c • The equation (9) then turns into the identity. 

In order to perform the perturbative expansion of the bound-state characte¬ 
ristics M and T around the unperturbed values we, as in Ref. [13], split the 
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full BS kernel V into two parts as V = Vc + V' and consider V' as a perturba¬ 
tion. Further, we introduce the projector onto the subspace orthogonal to the 
ground-space w.f. ijjc and the ’’pole subtracted” Coulombic Green’s function 


Q = l — N(M*)\ip c >< ipc\, G R (n = G c (n - i (10) 

It is easy to demonstrate that Eq.(9) can be rewritten in the following form 

< x\ =< x\N(M*)\ip c >< 4>c\[l + (AGo 1 - V')G R QY l (11) 

where AGq 1 = Gq 1 (P) — Gq 1 (P*). With the use of Eq.(ll) the following 
identity is easily obtained 

< 4> c \[l + (AGo 1 - V')G r Q}-\AG^ 1 - V')\ipc >= 0 (12) 

which is an exact relation and serves as a basic equation for performing the 
perturbative expansion for the bound-state energy. 

The equation (11) expresses the exact BS w.f. of the atom in terms of the 
unperturbed w.f. via the perturbative expansion in the perturbation potential 
V' . This potential consists of the following pieces: 

1. The purely strong part, which is isotopically invariant. This part survives 
when the electromagnetic interactions are ” turned off” in the Lagrangian. 

2. The part, containing the diagrams with the (finite) mass counterterms 
which are responsible for the m, r ± — electromagnetic mass difference. 

3. The part, containing the exchanges of one, two, ... virtual photons and 
an arbitrary number of strong interaction vertices. 

Note that the terms 1 and 2 are more important due to the following reasons. 
The first term includes strong interactions which are responsible for the decay 
of the pionium. The second term makes this decay kinematically allowed due 
to finite difference of charged and neutral pion masses. Consequently, it seems 
to be natural to consider together the pieces 1 and 2. We refer to the corre¬ 
sponding potential as V i2 - The T-matrix corresponding to the potential V i2 is 
defined by T V2 (P) = V ]2 (P) + Vi 2 (P)G 0 (P)Ti 2 (P). The rest of the potential 
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V' is referred as V 3 = V' — V\ 2 . In what follows we restrict ourselves to the first 
order in the fine structure constant a, i.e. consider the diagrams with only one 
virtual photon contained in V 3 . 

The ’’regular” part of the Coulombic Green’s function (10) can be split into 
two pieces according to GrQ = Go(M*)+5G. Here function 5G corresponds to 
the ladder of the exchanged Coulombic photons and thereby contains explicit 
powers of a. It is given by the following expression 


5G = i\Jw(p)w(q) 
$(p,q) = 167T m n a 


$(p, q) - S{p)S(q) 


8 d 


G 0 (M*,p)G 0 (M*,q) 


Up - q)- 


+ i R ( p,q) 


M * dM*. 

+ (m n a)~ 2 S(p)S(q)R(p,q) (13) 


where S/p) = 47rm T a0o(p 2 + p 2 ) \ R( p, q) = 25 - yj8/7rm n a[S( p) + S(q)]. 
The integral Ir(p, q) is given by 


^fl(p,q) = J — [p Up; P,q)--D 1 (0; P, q)], E* = ~m v a 2 (14) 

Returning to the basic equation (12) we expand it in the perturbative series 
considering V 3 and 5G as perturbations. Meanwhile we expand AG/ 1 in the 
Taylor series in SM = M — M* and substitute M = M* + A+ Ah/ 2 ) — 
i/2 Td) — i/2 Td) + (8M*) _1 Td) 2 + • • •. Restricting ourselves to the first order 
of the perturbative expansion we arrive at the following relation 

0 = -2 iM*SM- < ^c\Go\M*)G 0 (M)T 12 \^ c > (15) 

+ UU < MG«\M*)[G„(M) + Go(M)T 12 Go(M)] - gdy 1 Wc > 

- < V>cl(l + T l2 G 0 (M*))V 3 ( 1 + G„(M')T u Mc > 

dC (M*) 

- < ij>c\Go 1 (M*)[6M ^ J + G 0 (M*)T 12 G 0 (M*)\G^(M*)SGT 12 \7/;c > 

In the lowest order only the first two terms in Eq.(15) survive and in the 
second term Gq 1 (M*)Go(M) = 1 can be assumed. Then we obtain 
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If now we use the local approximation for T 12 , assuming that it does not depend 
on the relative momenta, we arrive at the well-known Deser-type formulae for 
the energy-level displacement and lifetime [5]. Note that on the mass shell 


Re(?T 12 ) ~ T(n + n —» 7 r + 7 r ), Im(?T 12 ) ~ w Am 7 r |T( 7 r + 7 r" 


ttV )! 2 


At the next step we assume that V 3 = 5G = 0 and evaluate the remaining 
integrals in Eq.(15), we arrive at the following result 


r( 2 ) 9 A 1 

—7-- = - 0.763a, Ei = — ma 2 (18) 

rd) 8 Ei ’ 4 v ; 

The first term of this expression called ” strong correction 1 ' was obtained in 
our previous paper [13]. However in difference with the present derivation in 
Ref. [13] we have used the Born approximation for the calculation of AE* 4 ), 
i.e. in Eq.(16) T 12 was substituted by Vi 2 . The last term comes from the 
relativistic normalization factor hi the instantaneous Coulombic 

potential (2). It arises since in the local approximation for the amplitude T 12 
the atom decay width is proportional to the quantity | / d^p/{2 'k) a, ^c{j > )\ 2 • For 
the particular choice of the potential (2) it is equal to 0 q(1 — 0.381 a) 2 /m n . 
Since this correction comes from the Coulombic w.f. of the atom, it does not 
depend on the parameters of the strong 7 T 7 t interaction, and for this reason it 
was neglected in Ref. [13]. Thus the name ’’strong” in Ref. [13] refers to the 
first-order corrections which survive in the limit when all terms containing 
explicit factor a (V3 and SG ) in the equation (15) as well as electromagnetic 
(relativistic) correction coming from the Coulombic w.f., are assumed to van¬ 
ish. 

The last term in Eq.(15), proportional to 5G, corresponds to the correction 
due to the exchange of the (infinite number) of Coulombic photons. With the 
use of the explicit expression for 5G (13) the integrals in this term can be 
easily evaluated, using again the local approximation for T 12 . These integrals 
are ultraviolet convergent, containing, however, an infrared enhancement alna 
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which stems from the infrared-singular one-photon exchange piece in Eq.( 6 ). 
Below we present the result of the calculations Q: 

_ ier(—2.694 + lira)—(19) 
m z Shi Ibnm 

Collecting all terms together and using Eqs.(16) for relating ImTi 2 to A E^\ 
we finally arrive at the following expression for the first-order correction to 
the 7 r + 7 r~ atom decay width 


r( 2 ) 

fw 


9 A Ed) 

8 Ei 
strong 


+ (—0.763a) + 

kernel normalization 


(1/2 +2.694-lna) 


A Ed) 

Ei 


Coulombic photon exchanges 


- ( tf * r m ) Re < i, c |(1 + T 12 G 0 (M*))U( 1 + G„(M*)T 12 )|*r > 


( 20 ) 


Let us now turn to the discussion of the last term in Eq.(20). The potential V 3 
present in this term contains the diagrams of two types: a) the diagrams where 
7 r + and 7 T~ lines are connected by photons, b) the diagrams where these lines 
are connected by strongly interacting particles (+ p, ...) as well. The diagrams 
of the first type correspond to the retardation correction [ 11 ], correction due 
to vacuum polarization [9], etc. The diagrams of the second type correspond to 
the radiative corrections [14]. In the Eq.(20) all these corrections are given in a 
closed form avoiding any difficulties connected with double counting problem. 
The kernel (1 + Ti 2 G 0 (M*))V 3(1 + Gq(M*)Ti 2 ) which appears in this term is 
constructed from the underlying Lagrangian with the use of the conventional 
Feynman diagrammatic technique. Note that the matrix element 


< i> c |(1 + T 12 G 0 (M*))V 3 (1 + Go(M*)T 12 )|^c > (21) 

can be written as < ^ 121 V 31 ?/q 2 > where ipi 2 = [1 + Gq{M*)Ti 2 ]^c stands for 
the solution of the BS equation with the kernel AV — V *- V 3 . The detailed 
analysis of the above mentioned corrections will be addressed in our future 
publications. 

1 This correction has been recently calculated by H.Jallouli and H.Sazdjian in 
the framework of the three-dimensional constraint equations (H.Sazdjian, private 
communication) 
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Apart of these corrections there are two types of corrections implicit in formu¬ 
lae (16). The first is clue to the finite difference of the masses of charged and 
neutral pions. In order to evaluate this correction one has to calculate the real 
and the imaginary parts of the ir + n~ —> 7r + 7r - transition matrix T 12 at the 
two-charged-pion threshold with the use of the low-energy chiral Lagrangian 
[4], in two cases: with and without the term responsible for the m,± - mj mass 
splitting. The second correction is caused by the dependence of the transition 
matrix T ]2 on the relative momenta of 7r + 7r - pair. This correction can be 
evaluated with the use of Eq.(16) provided the explicit expression of T i2 is 
known. 

In order to estimate the size of the first three terms in Eq.(20) we have used 
the following value of the singlet scattering length m n (2a,Q + ag) = 0.49 [19] 
corresponding to the value /S.E^/E\ = 0.24%. The first, second and third 
terms then contribute, respectively, —0.26%, —0.55% and +1.85%, to the de¬ 
cay width, and their total contribution amounts up to ~ 1%. The largest 
contribution comes from the cdna term in Eq.(20). Note that this logarithmic 
term is likely to cancel with the similar term coming from the last piece in 
Eq.(20) in analogy with the positronium case [15]. 
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